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Magnetohydrodynamic Jump Conditions for Oblique Relativistic Shocks 

with Gyrotropic Pressure 

Glen P. Double/ Matthew G. Baring,^ Frank C. Jones,^ and Donald C. Ellison^ 



ABSTRACT 

Shock jump conditions, i.e., the specification of the downstream parameters of the 
gas in terms of the upstream parameters, are obtained for steady-state, plane shocks 
with oblique magnetic fields and arbitrary flow speeds. This is done by combining the 
continuity of particle number flux and the electromagnetic boundary conditions at the 
shock with the magnetohydrodynamic conservation laws derived from the stress-energy 
tensor. For ultrarelativistic and nonrelativistic shocks, the jump conditions may be 
solved analytically. For mildly relativistic shocks, analytic solutions are obtained for 
isotropic pressure using an approximation for the adiabatic index that is valid in high 
sonic Mach number cases. Examples assuming isotropic pressure illustrate how the 
shock compression ratio depends on the shock speed and obliquity. In the more general 
case of gyrotropic pressure, the jump conditions cannot be solved analytically with- 
out additional assumptions, and the effects of gyrotropic pressure are investigated by 
parameterizing the distribution of pressure parallel and perpendicular to the magnetic 
field. Our numerical solutions reveal that relatively small departures from isotropy 
(e.g., ~ 20%) produce significant changes in the shock compression ratio, r, at all 
shock Lorentz factors, including ultrarelativistic ones, where an analytic solution with 
gyrotropic pressure is obtained. In particular, either dynamically important fields or 
significant pressure anisotropics can incur marked departures from the canonical gas 
dynamic value of r = 3 for a shocked ultrarelativistic flow and this may impact models 
of particle acceleration in gamma-ray bursts and other environments where relativis- 
tic shocks are inferred. The jump conditions presented apply directly to test-particle 
acceleration, and will facilitate future self-consistent numerical modeling of particle ac- 
celeration at oblique, relativistic shocks; such models include the modification of the 
fluid velocity profile due to the contribution of energetic particles to the momentum 
and energy fluxes. 
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1. INTRODUCTION 

Collisionless shocks are pervasive throughout space and are regularly associated with objects as 
diverse as stellar winds, supernova remnants, galactic and extra-galactic radio jets, and accretion 
onto compact objects. Relativistic shocks, where the shock speed is close to the speed of light, 
may be generated by the most energetic events; for example, pulsar winds, blastwaves in quasars 
and active galactic nuclei (e.g., Blandford & McKec 1977), and in gamma-ray bursts (e.g., Piran 
1999). They may naturally emerge as the evolved products of Poynting flux-driven or matter- 
dominated outflows in the vicinity of compact objects such as neutron stars and black holes. As 
relativistic shocks propagate through space, magnetic fields upstream from the shock, at even small 
angles with respect to the shock normal, are strongly modified by the Lorentz transformation to the 
downstream frame. The downstream magnetic fields arc both increased and tilted toward the plane 
of the shock and can have large angles with respect to the shock normal. Hence, most relativistic 
shocks can be expected to see highly oblique magnetic fields and oblique magnetohydrodynamic 
(MHD) jump conditions are required to describe them. 

Relativistic shock jump conditions have been presented in a variety of ways over the years. The 
standard technique for deriving the equations is to set the divergence of the stress-energy tensor 
equal to zero on a thin volume enclosing the shock plane and use Gauss's theorem to generate the 
jump conditions across the shock. For example, Taub (1948) developed the relativistic form of the 
Rankinc-Hugoniot relations, using the stress-energy tensor with velocity expressed in terms of the 
Maxwell-Boltzmann distribution function for a simple gas. de Hoffmann & Teller (1950) presented 
a relativistic MHD treatment of shocks in various orientations and a treatment of oblique shocks 
for the nonrelativistic case, eliminating the electric field by transforming to a frame where the flow 
velocity is parallel to the magnetic field vector (now called the de Hoffmann- Teller frame). Peacock 
(1981), following Landau k, Lifshitz (1959), presented jump conditions without electromagnetic 
fields, and Blandford & McKee (1976), also using the approach of Landau k. Lifshitz (1959) and 
Taub (1948), developed a concise set of jump conditions for a simple gas using scalar pressure. 
Webb, Zank, &; McKenzie (1987) provided a review of relativistic MHD shocks in ideal, perfectly 
conducting plasmas, and in particular the treatment by Lichnerowicz (1967, 1970), which used this 
approach to develop the relativistic analog of Cabannes' shock polar (Cabannes 1970), whose origins 
also lie in Landau &; Lifshitz (1959). Kirk k. Webb (1988) developed hydrodynamic equations using 
a pressure tensor, and Appl & Camenzind (1988) developed relativistic shock equations for MHD 
jets using scalar pressure and magnetic fields with components and (a parallel field with 
a twist). Ballard & Heavens (1991) derived MHD jump conditions using the stress-energy tensor 
with isotropic pressure and the Maxwell field tensor. By using a Lorentz transformation to the de 
Hoffman-Teller frame, they restricted shock speeds, uqj to uq/c < cosBb, where 0b is the angle 
between the shock normal and the upstream magnetic field and c is the speed of light; hence, this 
approach may only be used for mildly relativistic applications. 

All of these approaches assumed that particles encountered by the shock did not affect the shock 
structure, i.e., shocked particles were treated as test particles. Moreover, except for Kirk k Webb 
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(1988), they confined their analyses to cases of isotropic pressures, a restriction that is appropriate 
to thermal particles or very energetic ones subject to the diffusion approximation in the vicinity 
of non-rclativistic shocks. The assumption of pressure isotropy must be relaxed when considering 
the hydrodynamics of relativistic shocks, since their inherent nature imposes anisotropy on the ion 
and electron distributions: this is due to the difficulty particles have streaming against relativistic 
flows. The computed particle anisotropics for relativistic shocks are considerable in the shock layer 
(e.g., Bcdnarz & Ostrowski 1998; Kirk ct al. 2000), persisting up to arbitrarily high ultrarelativistic 
energies. The particles eventually relax to isotropy in the fluid frame far downstream, on length 
scales comparable to diffusive ones. However, only a small minority of accelerated particles achieve 
isotropy in the upstream fluid frame of relativistic shocks, due to the rapid convection to the 
downstream side of the flow discontinuity. These isotropized particles are present only when they 
manage to diffuse more than a diffusive mean free path. A, upstream of the shock; their contribution 
to the flow dynamics is therefore dominated by that of the anisotropic particles within a distance A 
of the shock. For non-linear particle acceleration, where the non-thermal ion or electron populations 
possess a sizable fraction of the total energy or momentum fluxes, calculating pressure anisotropy 
will be critical for determining the conservation of these fluxes through the shock transition (e.g., 
Ellison, Baring, &; Jones 1996). The jump conditions we develop here can serve as a guide to 
self-consistent solutions of non-linear relativistic shock acceleration problems (Ellison &; Double 
2002). 

Here, we extend previous work by deriving a set of fully relativistic MHD jump conditions 
with gyrotropic pressure and oblique magnetic flclds. We adopt the gyrotropic case as a specialized 
generalization because it (i) is exactly realized in plane-parallel shocks and is a good approximation 
for oblique shocks where the flow deflection is small, i.e., the fleld plays a passive role (generally high 
Alfvenic Mach number cases), and (ii) permits a comparatively simple expression of the Rankine- 
Hugoniot jump conditions. The sonic Mach number, Ms, used throughout our paper, refers to the 
shock speed compared to the last-mode magnetosonic wave speed as described in Kirk &; Duffy 
(1999), i.e., Ms = y^Spou^/ (5Po)) where po is the upstream mass density and Pq is the upstream 
pressure. The Alfven Mach number we refer to here is deflned as Ma = ^/ AttpqUq/ Bq (Bq is the 
upstream magnetic field), regardless of the shock Lorentz factor, 70 = [1 — (tto/c)^]~^/^; i.e., we use 
these definitions applicable to non-relativistic fiows as parameters for the depiction of our results 
at all 70. 

Our results are not restricted to the de Hoffmann- Teller frame and apply for arbitrary shock 
speeds and arbitrary shock obliquities. We solve these equations and determine the downstream 
state of the gas in terms of the upstream state first for the special case of isotropic pressure, and 
then, by parameterizing the ratio of pressures parallel and perpendicular to the magnetic field, for 
cases of gyrotropic pressure. 

A principal result of this analysis is that either dynamically important magnetic fields or sig- 
nificant pressure anisotropics produce marked departures from the canonical value (e.g., Blandford 
&; McKee 1976; Kirk &; Duffy 1999) of r = 3 for the shock compression ratio in an ultrarelativistic 
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fluid. The magnetic weakening of ultrarelativistic perpendicular shocks came to prominence in 
the work of Kennel & Coroniti (1984) on the interaction of the Crab pulsar's wind with its envi- 
ronment. The similarity of such consequences of fluid anisotropy and low Alfvcnic Mach number 
fields, which are also pervasive for trans-relativistic and non-relativistic shocks, has its origin in the 
similar nature of the plasma and electromagnetic contributions to the spatial components of the 
stress-energy tensor. This result may have important implications for the application of first-order 
Fermi shock acceleration theory to gamma-ray bursts and jets in active galaxies. 

In this work, we concentrate on using analytic methods for determining the fluid and elec- 
tromagnetic characteristics of the shock and do not explicitly include first-order Fermi particle 
acceleration. Future work will combine these results with Monte Carlo techniques (e.g., Ellison, 
Baring, & Jones 1996; Ellison &: Double 2002) that will allow the modeling of Fermi acceleration 
of particles, including the modification of the shock structure resulting from the backreaction of 
energetic particles on the upstream flow at all pertinent length scales. The jump conditions we 
present here, however, apply directly to test-particle Fermi acceleration in shocks with arbitrary 
speed and obliquity. 

2. DERIVATION OF MHD JUMP CONDITIONS 

2.1. Steady-State, Planar Shock 

Using a Cartesian coordinate system with the H-x-axis pointing towards the downstream direc- 
tion, we consider an infinite, steady-state, plane shock traveling to the left at a speed uq with its 
velocity vector parallel to the normal of the plane of the shock as shown in Figure 1. The upstream 
fluid consists of a thin, nonrelativistic plasma of protons and electrons in thermal equilibrium with 
TpQ = %.Q, where TpQ (Zio) is the unshocked proton (electron) temperature. A uniform magnetic 
field, Bq, makes an angle Q^o with respect to the x-axis as seen from the upstream plasma frame. 
We keep the field weak enough to insure high Alfven Mach numbers (i.e.. Ma > 2.5) and thus to in- 
sure that the magnetic turbulence responsible for scattering the particles is frozen into the plasma. 
The xyz coordinate system is oriented such that there are only two components of magnetic field, 
Bxo and B^o, in the upstream frame. The field will remain co-planar in the downstream frame and 
the downstream flow speed will be conflned to the x-z plane as well (e.g., Jones & Ellison 1987). 

In the shock frame, the upstream flow is in the -|-a;-direction and is described by the nor- 
malized four-velocity: /3g = ^0(1, f^Q ,0,0). The downstream (i.e., shocked) flow four-velocity is 
P2 ~ 72(1; /Sf' 0' /^l)' where the subscript (2) refers, here and elsewhere, to upstream (down- 
stream) quantities, /3 = u/c, and 7 = (1 — p'^)^^^'^ is the corresponding Lorentz factor associated 
with the magnitude of the flow three-velocity u. Note that these subscript conventions follow those 
used in Ellison, Berezhko, &; Baring (2000), where the subscript 1 is reserved for positions infinites- 
imally upstream of the subshock discontinuity, admitting the possibility of flow and field gradients 
upstream of the subshock. Here the use of subscript 1 is redundant. In addition, greek upper and 
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shock front 



Fig. 1. — Schematic diagram of a plane shock, viewed from the shock frame, showing obhque 
magnetic fields and a Gaussian volume over which the divergence of the stress-energy tensor is 
integrated. The angles and field strengths are measured in the local plasma frame and, in all of 
our examples, we take the upstream flow to be parallel to the shock normal. The y-axis is directed 
into the page. 
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lower indices refer to spacetime components 0-3, while roman indices refer to space components 
1-3. 

The set of equations connecting the upstream and downstream regions of a shock consist of the 
continuity of particle number flux (for conserved particles), momentum and energy flux conserva- 
tion, plus electromagnetic boundary conditions at the shock interface, and the equation of state. 
The various parameters that define the state of the plasma, such as pressure and magnetic field, 
are determined in the plasma frame and must be Lorentz transformed to the shock frame where 
the jump conditions apply. We assume that the electric field is zero in the local plasma rest frame. 
In general, the six jump conditions plus the equation of state cannot be solved analytically because 
the adiabatic index (i.e., the ratio of specific heats, whose value varies smoothly between the nonrel- 
ativistic and ultrarelativistic limits) is a function of the downstream plasma parameters, creating 
an inherently nonlinear problem (e.g., Ellison & Reynolds 1991). Even with the assumption of 
gyrotropic pressure, there are more unknowns than there are equations; however, with additional 
assumptions, approximate analytic solutions may be obtained. 



2.2. Transformation Properties of the Stress-Energy Tensor 

The stress-energy tensor, T'*'', describes the matter and electromagnetic momentum and energy 
content of a medium at any given point in space-time. Continuity across the shock is established by 

setting the appropriate divergence (or covariant derivative T^'^-iy) of the stress-energy tensor equal 
to zero, namely for directions locally normal to the shock. Following standard expositions such as 
Tolman (1934) and Weinberg (1972), the total stress -energy tensor, , is expressed as the sum 
of fluid and electromagnetic parts, i.e., 

-'tot ~ -'fluid -'em • \^) 

This is then integrated over a thin volume containing the shock plane as shown in Figure 1. Ap- 
plication of Gauss's theorem then yields the energy and momentum flux conditions across the 
plane of the shock by using T^^n-jy = 0. Accordingly, T^'^riy yields the conservation of energy flux, 
T^'^Uy yields the x-contribution to momentum flux conservation in the x-direction, T'^'^n^ yields 
the ^-contribution to momentum flux conservation in the x-direction, and riy = (0, 1, 0, 0) is the 
unit four-vector along the x-axis in the reference frame of the shock. The Einstein summation 
convention is adopted here and elsewhere in this paper. 

The components of the fluid and electromagnetic tensors are defined in the local plasma frame 
and are subsequently Lorentz-transformed to the shock frame where the flux conservation conditions 
apply, i.e.. 



-'tot.s — ^^a 



, (2) 



•■fluid,!! ^ -^EM,!! 

where the subscript n (s) refers to the plasma (shock) frame with the x-axis oriented normal to the 
shock in each case. Since the flow speeds in our model may have two space components, in the x- 
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and z-directions, the Lorentz transformation is 

/ 



A'' = 



\ 



7 


iPx 




7/3. 









1 





7/3. 







(3) 



where the P components and 7 are as previously defined in Section 2.1. The system defined in 
Figure 1 is invariant under translations in the y-direction, so that conservation laws along the 
y-axis are trivially satisfied. 



2.3. The Fluid Tensor and Equation of State 



The fiuid tensor will be constrained to the gyrotropic case in the local fiuid frame; i.e., pressure 
can have one value parallel to the magnetic field and a different value perpendicular to the magnetic 
field (with symmetry about the magnetic field vector) . This gives the diagonal stress-energy tensor 



2fluid,] 



/ e \ 

P\\ 

Pj_ 

\ P± J 



(4) 



where e is the total energy density, {P±) is the pressure parallel (perpendicular) to the magnetic 
field, and the subscript B refers to the magnetic axis in the plasma frame. We obtain the fluid 
tensor in the xyz plasma frame (labelled with subscript n referring to the shock normal), ?fiuid,n) 
with a rotation about the y-axis, 



uid,n 



uid,B 



7^- 



where 



7^ 



/ 1 \ 

cos ©B sin Bb 

10 

\ — sin @B cos Bb / 

The resulting tensor in the xyz plasma frame is 

/ e \ 

P|| cos^ Bb + Pl sin^ Bb (Pj_ - Py ) sin Bb cos Bb 

P± 

\0 (Pj_ - i^i ) sin Bb cos Bb P± cos^ Bb + Py sin^ Bb / 



(5) 



(6) 



^fluid,!! 



(7) 
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with the T*-' space components corresponding to the 3-dimensional pressure tensor presented by 
Elhson, Baring, & Jones (1996) for nonrelativistic shocks, once a typographical error is corrected. 

Renaming the components of the fluid tensor as 



^fluidm = 



/ 


e 


















p 





p 












p 

^yy 







V 





Pzx 





Pzz 


) 



(8) 



with Pxz = Pzx, we have the identities 



and. 



sin 8b 

' cos 6b ' 

cos Qb 
sin 6b ' 

Pxx 

+ 2P^, cot (263) 



Pi. — Pxx ~l~ Px 



(9) 
(10) 

(11) 



The T^^ component, e, is the total energy density in the rest or plasma frame. The other com- 
ponents, Pij, arc defined by Tolman (1934) as the "absolute stress" components in the proper 
frame. Pij is the pressure parallel to the i-axis exerted on a unit area normal to the j-axis. Hence, 
the diagonal components can be considered a pressure, but the off-diagonal components are shear 
stresses. The isotropic scalar pressure, P = rr(Pfluid,B)/3, is a Lorentz invariant. The fluid tensor, 
in general, changes its appearance significantly under a Lorentz transformation via the mixing of 
components: for example, it picks up momentum flux components in reference frames moving with 
respect to the proper frame (Tolman 1934). 

Using the above expressions, we can derive an adiabatic equation of state. Starting with the 
spatial portion of the fluid tensor: 



■Pfluid.B 





V 





p± 




\ 





(12) 



the total energy density can be written as 

_Tr(^fluid£) 



3(r - 1') + 



(13) 



where T'r(7'fluid,B) is the trace of the pressure tensor, F is the adiabatic index, and pc^ is the rest 
mass energy density. Using equations (9) and (10), 



^3(r-"lf^ = 3(r^_i) (^|| + 2^^^) = r^T + % (^^ot 6b - tan 63) 



(14) 
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In terms of the magnetic field, where tan 0b = Bz/B^, the adiabatic, gyrotropic equation of state 
becomes 



e = 



1 



r- 1 



Bz B-r 



(15) 



While r is well defined in the nonrelativistic and fully relativistic limits, in mildly relativistic shocks 
r depends on an unknown relation between e and P which we will address in a later section. 



2.4. The Electromagnetic Tensor 



The electric and magnetic field components of the general electromagnetic tensor in the plasma 
frame are given by Tolman (1934) as 



tEM,n 



47r 



E'^ + B^ 
2 


(E X B), 


(E X B)^ 


(E X B 


(E X B), 


Qxx 


Qxy 


Qxz 


(E X B)^ 


Qyx 


Qyy 


Qyz 


(E X B), 


Qzx 


Qzy 


Qzz 



where the Q's are the Maxwell stresses defined as 
Q 



e'^ + b'^ 



\ 



Ef - Bf and Q 



-iE^,Ej+BiBj) 



(16) 



(17) 



Note that here the suffix n denotes an axis orientation along the shock normal, as it does for the 
fluid tensor. Since the electric fields in the plasma frame are negligible, this simplifies to 



^5! 

2 



?EM,n - 



V 









2 




-BzBx 







B^ 
2 









—BxBz 



Bl - Bl 



\ 



(18) 



Observe that this electromagnetic contribution to the stress-energy tensor resembles the structure 
of that for the gyrotropic fluid in equation (8), i.e., the presence of laminar fields should mimic 
anisotropic pressures in terms of their effect on the flow dynamics. A noticeable difference, however, 
is that while the magnetic field can exhibit tension and can therefore generate negative diagonal 
components for T^*^, depending on the orientation of the field, the corresponding diagonal pressure 
components are always positive definite. 
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2.5. Flux Conservation Relations 



As discussed above, the energy and momentum conservation relations in the shock frame can 
be derived by applying T'^'^n^ = to equation (2), individually on the Lorentz-transformed fluid 
and electromagnetic tensors. The conservation of energy flux derives from 



^tot,s "-f — ^fluid,s ''^f + ^EM,s — ^ 

The fluid contribution to energy flux conservation is 



pfluid 
■ en 



7 [(27-1)/?^ + fP.z, 



while the electromagnetic contribution is 



TTiEM rpQu 



(19) 



(20) 



[{^-l)MlBl + {r^(5l + ^l)[i,Bl-{{2^-l)^l+^l}P,B,B,] . (21) 



The conservation of momentum flux derives from 



tot.s ' 



(22) 



The x-component of the transformed fluid tensor contributing to the conservation of momentum 
flux is 



Ft'" = Tfll„-id,s = 7'/3'(e + P^x) + Pxx - (7 - lf^{Pxx - Pzz) + 



px 



2(7-l)(7/3^+/3.')^P.. , 



and the z-component is 



Fpf" = T^luVs = 7'/9./?.(e + Pxx) - (7 - IM + iPD^iPxx - Pzz) + 



7 + 2(7-1)2^ 



(23) 



(24) 



The x-component of the electromagnetic contribution to the conservation of momentum flux is 
FT = r^M,s = ^PlB' + ^ m + Plf - (7 - {Bl - Bl) - 



^{^-l)i^(5l + (5l)^B,B,, (25) 
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and the z-component is 



FT = nVsn. = if^AB^ + ^(7 - - Pl)^{Bl - Bl) - 



^(7 - if^B^B, - ^B,B, . (26) 

Clearly, parallels between the various components of the fluid and electromagnetic stress tensors 
can be drawn by inspection of these forms for the fluxes. 

In all cases where the Alfven Mach number is greater than a few, the downstream flow velocity 
deviates only slightly from the shock normal direction so Pz Px- This allows a flrst-order 
approximation in f3z and the above equations become: 



Fen" ~ 7''Px{e + Pxx)+7i2j-l)PzPxz , 
pEM _ 7% 2 7(27 - 1) ,3 D D 

for the energy flux contributions, and 

F^x'" ~ ^^/3lie + Pxx) + Pxx + 2^{^-l)^Pxz , 



(27) 



g^PxB +g^{B,- B^) ^-——B^Bz , 



(28) 



and 



if;'<i ^ -f^PxPzie + Pxx) + iPxz - (7 - ^)%{Pxx - Pzz) 



(29) 



for the X— and z— components of momentum flux, respectively. As we show below, the approxi- 
mation, Pz *C Px, becomes progressively better as the shock Lorentz factor increases but, in fact, 
equations (27-29) provide an excellent approximation at all Lorentz factors for virtually the entire 
parameter regime we have considered in this paper. Unless Mach numbers less than a few and/or 
extreme anisotropics are considered, equations (27-29) yield solutions within one part in 10"^ to 
those obtained with equations (20-26) and are much easier to solve. 



2.6. Jump Conditions 



The jump conditions consist of the energy and momentum flux conservation relations, the 
particle flux continuity, and the boundary conditions on the magnetic fleld. The conservation of 
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particle number fiux^ is 



n 2 



0, 



where the brackets provide an abbreviation for 



(30) 



(31) 



This jump condition, as well as the ones that follow, are written in the shock frame and, as always, 
the subscript (2) refers to upstream (downstream) quantities. The remaining jump conditions 
are: 

1 2 

= 0, 



pfluid I pEM 

en en 



F- 



px 



EM 
px 



F- 



fluid 



pz 



+ F, 



EM 
pz 



J 

- 2 

- 

2 

J 



, 



= . 



(32) 



Adding the steady-state conditions on the magnetic field, namely that V • B = 0: 



and also that V x E = 0: 



Bx =0 , 



-fiP.Bx - PxB, 



= 0, 



(33) 



(34) 



completes the set of six jump conditions. In the limit of (7 — 1) ^ 1 and for isotropic pressures, 
the above expressions reproduce the standard continuity conditions at non-relativistic shocks (e.g., 
see p. 117 of Boyd & Sanderson 1969). 

At this point there are eight unknown downstream quantities (Px2, Pz2, Bx2, Bz2, Pxx2, Pxz2, 
62, and 712) and only six equations. If isotropic pressure is assumed, Pxx = Pzz = P and Pxz = 0, 
leaving six equations and seven unknowns. To obtain a closed set of equations for isotropic pressure, 
an assumed equation of state (e.g., equation 15) is added to the analysis. Successive elimination 
of variables then generally leads to a 7th-order equation in the compression ratio r = f3xo/l3x2 
with lengthy algebraic expressions for its coefficients (e.g., Webb, Zank, & McKenzie 1987; Appl 
&L Camenzind 1988). The forms of the coefficients depend on the assumed equation of state and 
do not simplify easily. Here we perform some of the simpler algebraic eliminations and then use a 
Newton-Raphson technique to iteratively solve two simultaneous equations. 

In Section 3.1 below, we derive an approximate expression for the downstream adiabatic index, 
r2, for cold upstream plasmas. For anisotropic cases, further microphysical information is required. 



*We assume there is no pair creation nor annihilation. 
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which is generally only accessible using computer simulations; in the gyrotropic approximation, we 
parameterize pressure anisotropy in Section 3.2 to provide insight into global characteristics of the 
jump conditions. 



3. RESULTS 

3.1. Isotropic Pressure, High Sonic Mach Number Cases 

Oblique shock jump conditions cannot, in general, be solved analytically even for isotropic 
pressure because the downstream adiabatic index, depends on the total downstream energy 
density and the components of the pressure tensor (or scalar pressure), which are not known before 
the solution is obtained. The problem is inherently nonlinear except in the nonrelativistic and 
ultra-relativistic limits where r2 = 5/3 and 4/3, respectively. Furthermore, the gyrotropic pressure 
components are determined by the physics of the model and do not easily lend themselves to analytic 
interpretation, although Kirk & Webb (1988) provided equations based on a power-law distribution 
in momentum for the pressure tensor components in the special case of a parallel relativistic shock 
with test particle first-order Fermi shock acceleration. 

An excellent approximation can be obtained in the absence of efficient particle acceleration if 
^^th ^ ^0) where Vth is the thermal speed of the unshockcd plasma. In this case, which corresponds 
to high sonic Mach numbers, upon scattering in the downstream frame all particles have 

7rel^(l-^.'el)~'/', (35) 

where 

/^rel = — = ^ (36) 

C i — P0P2 

is the relative /3 between the converging plasma frames. Prom kinetic theory, the isotropic pressure 
is 

Ti 

P = -<p.v>, (37) 

where n is the particle number density, and p and v are the particle momentum and velocity, 
respectively. Then, with our approximation for particle velocity and using the isotropic version of 
equation (13), 

V2 = ^ + 1 = (V^)^^^'^^^l + 1 = 7rel/?rel ^ (3g) 
^ e - (7rel - l)mc2 3(7rel - 1) 

or, 

= . (39) 

37rel 

This essentially kinematic approximation, which is operable only if diffusive transport of particles 
from downstream to upstream contributes insignificantly to the momentum and energy fluxes, per- 
mits a direct numerical solution for isotropic pressure, arbitrary obliquity (as long as the upstream 
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M, > 100 



Q =5° 



> 100 



1 — I I 1 1 1 I I 1 1 1 — I I I iiiij 1 — I I I 

M. = 60 (6) 




0.01 0.1 1 10 100 

^oyo 



0.1 1 10 100 



Fig. 2. — Compression ratio, r, @b2, and Q^2 versus /?o7o- The three panels on the left have 
Ma ^ 100 and the sonic Mach numbers as shown. The three right-hand panels have Ms > 100 
with the Alfven Mach numbers as shown. In all cases, 9bo = 5°. Note that in the top right-hand 
panel, the Ma = 6 and 60 curves are nearly identical. 
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Alfvenic Mach number, Ma, is high), and arbitrary flow speed (see, for example. Kirk 1988; Gallant 

2002, for alternative forms for r2). Note that equation (39) provides an upper limit to the adiabatic 
index because any particles accelerated by the shock would tend to raise the average Lorentz factor 
and cause the adiabatic index to decrease. 

In Figures 2 and 3 we show results for the compression ratio, r, and the downstream angles, 
Ob2 and Ou2, as a function of Po'jo, for two extreme upstream magnetic field angles, 6bO) and 
various sonic and Alfven Mach numbers. The magnitude of the downstream field is 

S2 = [5^0 + 7o'(r2- 1)5^0] V2. (40) 

There are a number of important characteristics of these results. First, the jump conditions 
map smoothly from fully nonrelativistic to ultrarelativistic shock speeds and obtain the canonical 
values for the compression ratio r = Pxo/Px2 = 4 for high Mach number, nonrelativistic shocks, and 
r = 3 for high Mach number ultrarelativistic shocks. For /3o7o ^ 1 and regardless of the obliquity, 
a low Ms results in a weaker shock with smaller r, as expected (similar behavior is exhibited in 
Figure 1 of Appl & Camenzind 1988). For /3o7o ^ 10, the sonic Mach number Ms has little influence 
on the results until it becomes very low, i.e. 7oMg 1. Since our definition of Ms is inherently 
non-relativistic, perceptible changes to the fluid dynamics arise only when the pressure Pq becomes 
comparable to the relativistic ram pressure 7o/3oPoc^; this domain is exhibited in the parallel fluid 
shock jump conditions explored by Taub (1948). 

Figure 2 with 0bo = 5° and Figure 3 with 9bo = 85° illustrate an important characteristic of 
relativistic shocks. For all upstream field obliquities other than 6bo = 0, the downstream magnetic 

field angle shifts towards Bb2 = 90° as the shock Lorentz factor increases, indicating the importance 
of addressing oblique fields when treating acceleration at highly relativistic shocks. The transition 
criterion is directly obtainable from equation (34) (since 0^2 is, in general, very small). One 
quickly arrives at /3o7otanOBo ^ 1 being the necessary condition to render tan6B2 > Ij so that 
0B2 90°. Furthermore, in the left hand panels of Figures 2 and 3, the angle the downstream 
flow makes with the shock normal, 0u2, is small at all 7o/3o (note the logarithmic scale for 6u2), 
consistent with the assumption that Pz ^ Px- This is a consequence of the passive magnetic fleld 
corresponding to Ma ^ 1. Very different behavior is exhibited when the fleld becomes dynamically 
important, as is evident in the upper right-hand panels of Figures 2 and 3: a low A/a (i.e., high 
-Bo) can drastically lower r, even at ultrarelativistic speeds (see Kennel &: Coroniti 1984). 

In Figure 4 we show r2 for the two extreme sonic Mach number cases from Figure 2. Our 
approximation for r2 depends on uq and r, but turns out to be quite insensitive to r, at least in 
the range above r ^ 2.7. The r2 generated by equation (39) is essentially the same as the low 
temperature solutions presented in Figure 2b of Heavens & Drury (1988) for an e+e~ plasma. As 
noted above, if Fermi acceleration is permitted to occur, r2 will approach 4/3 at lower /3o7o due to 
the contribution of energetic particles. 

The variation of r with shock speed we show here closely resembles the low sonic Mach number 
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Fig. 3. — Same as in Figure 2 except all examples here have Oro = 85°. Note the drop in r, even 
at 7o ^ 1, in the upper right-hand panel where anisotropic magnetic stresses become important. 
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Fig. 4. — The downstream ratio of specific heats, T2, versus /3o7o for the extreme Afs cases shown 
in the left-hand panels of Figure 2. Very similar curves are obtained for all of the examples shown 
in Figures 2 and 3. 
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solutions depicted in Figure 4 of Fujimura & Kennel (1979), who considered jump conditions in 
parallel (i.e., Ma oo) trans-relativistic shocks using an isotropic Jiittner-Synge equation of state 
(Syngc 1957). The principal effect of upstream heating is to lower the compression ratio and 
weaken the shock when the ratio of the particle pressure to the rest mass energy exceeds the square 
of the shock four-velocity. The compression ratio clearly becomes insensitive to the plasma heating 
in ultrarelativistic, parallel shocks since Fq is always 4/3 and the solution is /3o/32 = 1/3 (e.g., 
Blandford &; McKee 1976). An array of possible jump conditions is admitted when an extension 
to multi-component plasmas is explored, such as in Peacock (1981), Kirk (1988), Heavens & Drury 
(1988) and Ballard & Heavens (1991), where the thermal interplay of ions and electrons on the shock 
dynamics in the trans-relativistic regime can be encapsulated using the two adiabatic shock index 
parameters Fq and F2. Such a parameterization implies that a variety of equations of state can be 
accommodated within the formalism presented here. Notwithstanding, when particle acceleration 
is considered, thermal equations of state become inappropriate and simulation results become more 
important (Ellison & Reynolds 1991). 
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Fig. 5. — Compression ratio versus Alfvcn Mach number for various upstream magnetic obliquities, 
0BO- The values of r are obtained in the ultrarelativistic limit with a = 1, 
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Figure 5 shows how r varies as a function of Ma in the ultrarelativistic limit for various values of 
0BO cind with isotropic pressure. The compression ratio depends strongly on the upstream obliquity 
and can drop well below 3 for low Ma- The lower right-hand panels of Figures 2 and 3 also show 
that a low Ma produces magnetic stresses that cause the downstream flow to deflect from the 
shock normal direction, causing to vary inversely as Ma- Despite the fact that 0u2 ~ 20° for 
Ms = 100 and Ma = 2.6, the approximate equations (27)-(29) give essentially identical results as 
the complete equations (20-26). The impact of either Mach number on 6b2 is relatively small. 

The most important consequence of a large magnetic field energy density is that it lowers the 
compression ratio in oblique shocks (the field is dynamically passive in parallel ones with 0bo = 0) , 
even when 70 S> 1. The effect is contained in Eq. (4.11) of Kennel &; Coroniti (1984), which specifies 
the jump condition or downstream flow four- velocity for an ultrarelativistic perpendicular MHD 
shock. Algebraically simplifying their formula, and expressing it in terms of the three-velocity 
compression ratio r and Alfvenic Mach number Ma used here, leads to the form (for (3x0 = /3o ~ 1) 



where the ratio a of the magnetic plus electric energy flux to the particle energy flux that is used 
by Kennel &; Coroniti (1984) is given by a « 1/M| = Bq / {An pqUq) . This result apphes to 70 » 1 
regimes, and is reproduced in the numerical results depicted in the top right hand panel of Fig 3 
and also the 89° curve of Fig 5. 

The weakening of nonrelativistic shocks at low Alfven Mach numbers is widely understood. 
Such an effect is suggested for shocks with f3xo ^ 0.8 in Ballard & Heavens (1991) and is also 
somewhat apparent for P^o ^ 0.97 in Appl &: Camenzind (1988). Kirk &; Duffy (1999) show r 
as a function of Ma through the trans-relativistic regime for 6bo = 45°. The mildly relativistic 
regime was appropriate for shocked jets in active galactic nuclei, the main application of relativistic 
MHD shock analyses over a decade ago. Here, it is evident that this lowering of r by dynamically 
important magnetic fields persists up to arbitrarily high 70, a result obtained by Kennel & Coroniti 
(1984), who applied relativistic MHD to the consideration of perpendicular pulsar termination 
shocks. The more recent association of ultrarelativistic shocks with cosmic gamma-ray bursts 
further motivates the extension to the 70 > 100 regime. Accordingly, this magnetic weakening 
of the shock has profound implications for the interpretation of gamma-ray burst spectra and 
associated emission mechanisms. 

The origin of this effect can be attributed to the anisotropic and intrinsically relativistic nature 
of the field structure. In the ultrarelativistic regime, the equation of state of a quasi-isotropic, 
turbulent field structure replicates that of the familiar ultrarelativistic gas with T = 4/3. However, 
if the field is laminar and oblique, the stress-energy tensor in equation (18) exhibits anisotropic 
stresses. This anisotropy influences the flux conservation relations for dynamically important fields 
when Bzo 7^ 0; moreover it should mimic to some extent the behavior anticipated from anisotropies 
due to the gas contribution. Such a parallel obviously emerges from results presented in the next 
Section. 
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Fig. 6. — The top three panels show the compression ratio, r, versus /3o7o for shoeks with varying 
downstream anisotropy, a = P±/P^^, and obhquity, Gbo- In ah panels, the solid curves have a = 1, 
the dashed curves have a = 0.8, and the dotted curves have a = 1.2 downstream. The upstream 
plasma is taken to be isotropic (i.e., a = 1) in all cases, but the upstream a has virtually no 
influence on the solutions for the parameter regime we consider. The anisotropy has little effect on 
Qb2, but does cause modifications in 0u2 at mildly relativistic shock speeds, as indicated in the 
bottom panel which shows the angle the downstream flow makes with the shock normal for the 
particular case 6bo = 10°. These are high Mach number examples (Ms = Ma > 100); the effect 
on r at lower Mach numbers is shown in Figure 7. 
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3.2. Gyrotropic Pressure 

For gyrotropic pressure, an additional constraint is needed to close the system of equations and 
obtain an analytical solution; we impose this via 

P± = aP\\ , (42) 

where a is an arbitrary parameter and a = 1 corresponds to isotropic pressure. Equation (42) 
allows us to illustrate the effects of anisotropic pressure, but is not suggested as a model for 
specific acceleration scenarios. In relativistic shocks where the accelerated population contributes 
significantly to the total dynamical pressure, the value of a should deviate significantly from unity 
in both the upstream and downstream regions. Using equation (7), equation (42) then yields 

Pxx = P\\ (cos^ Bb + a sin^ Gb) , (43) 

Pxz = P\\ (a - 1) sin Gb cos Gb , (44) 

and using the fluxes given in equations (20)-(26) or (27)-(29), we have a closed set of equations 
for the jump conditions for shocks with gyrotropic pressure, arbitrary obliquity, and arbitrary flow 
speed. Results for various Gbo's and a's are shown in Figure 6 (these results all have Ms = Ma = 

100 but they remain unchanged for larger Mach numbers). The solid curves have a = I, the dashed 
curves have a = 0.8, and the dotted curves have a = 1.2. In all cases, we have taken the pressure 
in the unshocked gas to be isotropic and a ^ 1 is only applied downstream. 

While our solutions can allow for anisotropic upstream pressure, an upstream a 7^ 1 generally 
produces insignificant changes to our results for the parameter regime discussed here. Notable 
exceptions arise when the sonic Mach number is very low, i.e. 7oM| ~ 1. We observe that angular 
distributions at relativistic shocks generated by Monte Carlo simulations (e.g., Bednarz &z Ostrowski 
1998; Ellison & Double 2002) and semi-analytic convection-diffusion equation solutions (Kirk et al. 
2000) indicate that for plane-parallel scenarios with Gbo = the accelerated particles are dominated 
by parallel pressure upstream (a < 1) and perpendicular pressure downstream (a > 1) when near 
the shock. This would suggest a possible weakening of the shock if the non-thermal population 
were to contribute significantly to the dynamics, thereby rendering such a contribution less likely. 
The picture may be much different for oblique and perpendicular shocks. 

The eflFects of anisotropic pressure on the compression ratio depend strongly on 70 and Gbo- 
When Gbo is small, the downstream angle Gb2 is also small at nonrelativistic and mildly relativistic 
shock speeds (top panel of Figure 6). Therefore, at these speeds 

Pxz ~ and Pxx ~ P\\ for any 

a. Since = P±/a, the fraction of downstream pressure in P^x is inversely proportional to a 
and since Pxx largely determines r, the compression ratio is less than the isotropic value, riso? fc*r 
a = 0.8 and greater than rjso for a = 1.2, as shown in the figure. In contrast, as 70 — >■ 00, Qb2 
approaches 90° for any Gbo > (see Figure 2) and Pxx ~ P± with Pxz again approximately equal 
to zero. Now, the fraction of pressure in Pxx will be proportional to o and r <! t\so for o; > 1 
and r > Viso for a < 1. The transition where r crosses riso occurs at slower shock speeds as Gbo 
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Fig. 7. — Compression ratio versus anisotropy for ultrarelativistic shocks. The top panel shows the 
effects of ©BO) while the bottom panel shows the effects of Mach number keeping Obo = 89°. Note 
that, depending on a = Pi/Py, the compression ratio r can be either larger or smaller than 3. 
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increases (see the Bbo = 10° panel of Figure 6) until 6bo is large enough (6bo = 60° panel) so no 

transition occurs. While the examples in Figure 6 all have Ms = -Ma = 100, the transition /3o7o 
is independent of Mach number. The dependence of 0b2 on a is relatively small for the examples 
shown in Figure 6 but Ou2 can change significantly with anisotropy, as shown in the bottom panel 
for the 6b2 = 20° example. 

In Figure 7 we show r as a function of a for various Mach numbers and Obo's, all in the 
ultrarelativistic limit. The top panel shows that r is relatively insensitive to the upstream magnetic 
field angle, with the lower panel showing a somewhat greater sensitivity to Mach number. The most 
important aspect of these plots is the fact that r can be either higher or lower than the canonical 
value of 3 if anisotropic pressure is important. This contrasts with the effects of a magnetic field 
(with isotropic pressure), where low Ma's in the ultrarelativistic limit gave r < 3 only (Figures 3 
and 5). 



3.2.1. Analytic solution for 70 S> 1 

Using our flux conservation relations in the limit 70 ^ 1, and assuming that the upstream 
pressure is isotropic and the downstream adiabatic index r2 = 4/3, we obtain an analytic solution 
for r valid for any oblique angle and any Ms or Ma greater than one: 

(a + l)bor^ + 2a{wo + boy - [2(3a + l)wo + (7a + 3)bo] r + 2(2a + l){wo + 60) = , (45) 

where wq = eo + Pq is the enthalpy density and 60 = Bq^/Att. Only the z-component of B appears 
since 3^2 = B^o and this component drops out. Eliminating the r = 1 root gives 

(a + l)6or^ + [2awo + (3a + l)6o] r - 2(2a + \){wo + bo) = , (46) 

or, 

\[2aq + (3a + 1)]^ + 8(a + l)(2a + l){q + l)}^^^ - 2aq - (3a + 1) 

^ = - 2W+T) ' 

where q = wo/bo- Solutions in terms of Ma and Gbo can be found using the relation 60 = 
Pqc^ sin^ Qbo/MJ^. For isotropic pressure, i.e., a = 1, 

(| + l)' + 3(g + l)j'^'-(| + l) . (48) 

Equations (48) and (47) reproduce the results obtained with the exact equations shown in Figures 5 

and 7 (bottom panel, solid curve) to a high degree of accuracy. In the limit of Afs ^ 1 and 
©BO ~^ 90°, equation (48) is equivalent to equation (41) from Kennel & Coroniti (1984). Once r 
is determined, the downstream z-component of B can be found from equation (40). 
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4. CONCLUSIONS 

We have derived shock jump conditions for arbitrary shock speeds and obhquities. When 
combined with a simple approximation for the ratio of specific heats (equation 39) , these equations 
specify the downstream conditions in terms of upstream parameters for isotropic pressure. For the 
case of gyrotropic pressure, an additional arbitrary parameter (equation 42) is required to close the 
set of equations and we have presented a number of solutions where the downstream pressure is 
gyrotropic. The exact equations for energy and momentum conservation are fairly complicated, but 
we have presented simpler, approximate results (equations 27-29) which are extremely accurate in a 
wide parameter regime (i.e., Ms > 5; M\ > 5; 0.8 < a < 1.5) for any shock speed or obliquity. To 
our knowledge, this is the first presentation of oblique shock solutions with gyrotropic pressure that 
continuously span the domains of nonrelativistic, trans-relativistic, and ultrarelativistic shocks of 
arbitrary obliquity. In addition, we have presented an analytic solution for the shock compression 
ratio, r, in the limit of ultrarelativistic shock speeds valid for oblique shocks of any Mg or Ma with 
gyrotropic pressure (equation 47). 

The results presented here assume that no first-order Fermi acceleration occurs, but they ap- 
ply directly to test-particle acceleration where the energy density in accelerated particles is small. 
They also constitute an important ingredient in more complex models of nonlinear particle accel- 
eration. In the test-particle case, the compression ratio r is altered by large magnetic fields and/or 
anisotropic pressures, even at ultrarelativistic speeds. The observation of marked departures from 
the canonical value of r = 3 for the compression ratio of an ultrarelativistic shock, for either dy- 
namically important fields or significant pressure anisotropics, is a major result of this paper. In 
the case of magnetic field influences on the dynamics of relativistic shocks, our results extend the 
conclusions of Kennel & Coroniti (1984) to all shock obliquities. The similarity of such consequences 
of fluid anisotropy and low Alfvenic Mach numbers is a consequence of the similar nature of the 
plasma and electromagnetic contributions to the spatial components of the stress-energy tensor. 

The importance of this result is obvious, since changes in r map directly to changes in the 
power-law index of the accelerated spectrum. This power law is the most important characteristic 
of test-particle Fermi acceleration, and is usually associated with trans-relativistic internal shocks 
in gamma-ray burst (GRB) models (e.g., Rees &: Meszaros 1992; Piran 1999). In addition, the outer 
blast wave is an ultrarelativistic shock during most of its active phase, sweeping up and accelerating 
interstellar material. This shock, believed to produce long-lasting afterglows, eventually transitions 
to a non-relativistic phase. Our results are directly applicable to test-particle acceleration models 
of both the internal and external shocks in GRBs, as well as to shocks believed to exist in jets 
in active galactic nuclei. In future work, we will apply these jump conditions to nonlinear shock 
acceleration models where the accelerated particle population can modify the shock structure. 
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